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Abstract 

We investigate as a member of the Ising universality class the Next-Nearest 
Neighbour Ising model without external field on a simple cubic lattice by 
using the Monte Carlo Metropolis Algorithm. The Binder cumulant and the 
susceptibility ratio, which should be universal quantities at the critical point, 
were shown to vary for small negative next-nearest neighbour interactions. 

keywords: next-nearest neighbour Ising model, Binder cumulant, fourth order cu- 
mulant. 

1 Introduction 

The great variety of critical behaviour is reduced by dividing all systems into a small 
number of universality classes. These classes are characterized by global properties 
like dimensionality and the number of components of the order parameter and have 
the same set of critical exponents and scaling functions, which are independent of 
microscopic details [l]. Here the 3D Ising universality class is considered by studying 
the Binder cumulant [2j of the Next-Nearest Neighbour Ising model (NNN model). 
If all universal quantities of one class are known, the asymptotic critical behaviour is 
believed to be also known completely under the condition, that only two nonuniversal 
amplitudes Kt and are given. This particular feature is called two-scale factor 
universality. The free-energy-density /, from which we can derive thermodynamical 
quantities and which is measured in units of ks, can be expressed by 0: 

f{t,H-L) = Ut,H;L) + Us{t;L), 

where t = 1 — T/T^ and H is the ordering field. The singular part fs is the part 
which yields thermodynamical singularities in the L ^ oo limit, fns denotes the 
non-singular part of /, which can be chosen without field dependence. For the 
singular part fs we can write for the Ising model 

fsjt, H- L) = L-''Y{KttL^I\ K^HL"'') 

Y{x, y) is called the scaling function and is universal using system-dependent factors 
Kh and Kt. 

Chen and Dohm ^ recently predicted some deviations from this universality, from 
a (p^ theory. 

In the following section we study the universal value of the Binder cumulant for 



the NNN model for different next-nearest neighbour interaction Jnnn by simula- 
tions. According to Privman, Hohenberg and Aharony |3j, the Binder cumulant 
scales according to 

U^it,L)^G{KttL'/'') , 

where G{x) = [{d^Y/dy^)/{d'^Y/dyy]y=o- For t = (means T = T^), f/4(0,L) 
approaches with growing L a universal constant Ul{0,oo), U^{0,L) which we are 
going to simulate, shows only small deviations from the universal value ^74(0, cxd) 
[S]. Consequently, it is possible to give a realistic approximation of the final value 
f/4 (0, 00) simulating finite lattices. 

2 Applied NNN Models 

2.1 Standard Next-Nearest Neighbour Model 

To show the variation of within the 3D Ising universality class, we choose the 
Next-Nearest Neighbour model (NNN Model) without external field and a simple 
cubic lattice. Each lattice site has 6 next neighbours in the distance a and 12 next- 
nearest neighbours in the distance ^/2- a, if a is defined as the lattice constant. The 
exchange force Jnnn of the next-nearest neighbours can now be chosen to be differ- 
ent from the exchange force Jnn of the 6 nearest neighbours as we also would also 
intuitively expect, because their distance to the considered site is also not identical. 
A ferromagnetic phase transition is also existent for an antiferromagnetic exchange 
force Jnnn- This kind of exchange, ferromagnetic for nearest neighbours and anti- 
ferromagnetic for next-nearest neighbours, will be treated for our universality test 
off/4*. 

2.2 Anisotropic Next-Nearest Neighbour Model 

Furthermore we studied the universality of in the zero field NNN Ising Model 
with an isotropic nearest-neighbour (NN) coupling J > and an anisotropic next- 
nearest neighbour (ANNN) coupling Jnnn < 0. We refer to the considered model 
of Chen and Dohm in their publication "Nonuniversal finite-size scaling in anistropic 
systems" [4]. The anisotropic NNN Ising Model is established by considering only 6 
of the 12 next-nearest neighbours being effective for NNN interaction Jnnn < and 
the other 6 NNN having no NNN interaction. The effective next-nearest neighbours 
have the positions ±(1, 0, 0), ±(1, 0, 1) and ±(0, 1, 1) of a simple cubic lattice. 

3 Simulation Results 

Simulating the NNN model Jnnn can be varied, while it is chosen Jnn = 1- For val- 
ues of Jnnn higher than a critical value JNNN,crit ~ —0.5 we observe ferromagnetic 
phase transitions, while below Jnnn^ctu no ferromagnetic phase transition occurs. 
At Jatatat = we obtain the standard 3-dimensional (3D) Ising ferromagnet. Here 
the region JNNN,crit < Jnnn < and thus a small antiferromagnetic interaction 
Jnnn is considered. 



Simulating now the "fixed point value" at the critical point of the Binder cumulant 
= U^iTc), we evaluate U^{Jnnn) to check its universality within the 3D Ising 
universality class. For an Ising Model in zero field, the Binder cumulant simplifies 
to 

U. = l~^, (1) 

where M symbols the magnetization per spin. For the determination of each f/4 
was plotted f/4 as a function of temperature T for three different lattice lengths L 
(= 10, 20, 40). To get more precise results, several curves of Ul{T) for different 
lattice sizes with different seeds were averaged. In figure 1 is shown this method for 
two different Jnnn- The errors were calculated by averaging the values of the three 
different crossing points of L = (10,20), (10,40), (20,40). 
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Figure 1: Localization of the fixed point {Ul,Tc) for different Jnnn by quantifying 
the crossing point of U^iT) for lattice sizes L = 10 (dotted), 20 (dashed) and 40 
(solid). 



3.1 Results for the Standard NNN Ising Model 

In figure El and table 1 we can see all simulated results as a function U^{Jnnn) and 
these results in comparison with theoretical results of Chen and Dohm |lj plotted 
as a function 1 — U^{w)/U{0) with w = Jnnn/{Jnn + '^Jnnn)- Chen and Dohm 
considered an anisotropic NNN Ising model, which we will discuss later. It clearly 
can be seen in figure and, that is not a fixed value for different Jnnn, particularly 
in the interval of Jnnn = [-0.25:-0.22]. This implies in this region the absence of 
universality for the 3D Standard NNN Ising model. 

Great deviations of the resulting f/4 for L — > 00 can be excluded for Jatattv = (H]. 
Observing lattice length up to L = 100 for Jnnn 7^ we observed no significant 
change of our simulation results. 

We can conclude, that in the region of small negative values for Jnnn the 3D Next- 
Nearest Neighbour Ising model shows a different universal behaviour than the 3D 
Nearest Neighbour Ising model. 
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Table 1: Simulation results of the universal quantity f/4 and its statistical errors for 
different next-nearest neighbour interaction Jnnn- 
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Figure 2: Standard 3D NNN model. Left side: Simulation results of the fixed 
point f/4 as a function of Jnnn', Right side: Results of the figure on the left (■) 
in comparison with theoretical results (for a different neighbourhood) of Chen and 
Dohm |1| (dashed line). 



3.2 Results for the ANNN Ising Model 

Figure inland table 2 show the results for the ANNN Ising Model analogue to section 
13.11 The simulation results show a significant non-universality of the Binder cumu- 
lant. As compared to the analytic results of Chen and Dohm ^ the analytical curve 
have a different tendency and is not within the errorbars of the simulation results. 
For w —0.5, f/4 increases whereas Chen and Dohm obtained a decrease. Thus, 
the theoretical results cannot be confirmed fully. 
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Figure 3: Left side: Simulation results (+) for the ANNN model of 
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side: Simulation results ofl—Ul{w)/U^{0) in comparison with the theoretical results 
(for this neighbourhood) of Chen and Dohm |ij (solid line). 
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Table 2: Simulation results for the ANNN Model for Ul{w). 



3.3 Universality Test of the Susceptibility for the ANNN 
Model 

Simulating the susceptibility x = (N/kBT) ((M^) - {\M\f) for the ANNN Ising 
Model we use the same range of NNN interaction Jnnn- We consider two temper- 
atures near the critical temperature Tc, T_ = Tc — 3%Tc and T+ = Tc + 3%Tc. The 
results show explicitly the absence of universality. Figure |31 shows the ratio of the 
susceptibility at T+ and T_ as a function of w. A significant deviation of nearly 100% 
exists here between w = and w = —0.45. Varying below Tc the lattice size L be- 
tween 10 and 317 and the number of Monte Carlo steps MCS between 50 and 5 ■ 10^ 
we can exclude finite-size effects with L = 40 and MCS = 5 ■ 10^. Table El shows 
all results of x{T-) and x(T+)/x(T_) for different Jnnn- (We looked at the 

fiuctuations in the absolute value of the magnetization < > — < \M\ not of 
the magnetization. Thus our values above are lower by a factor 1—2/tt = 0.3634 
than the usual < >, assuming Gaussian fiuctuations |2|.) 
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Figure 4: Simulation results of the ratio of susceptibility x{T+)/x(T^) as a function 
of Jnnn- 
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Table 3: Simulation results of the susceptibility at T_ = — 3%Tc and T+ = 
Tc ± 3%Tc and of the ratio x{T+)/x{T~) as a function of Jnnn- 



4 Conclusions 



This work examined the universality of the Binder cumulant for the 3D NNN Ising 
model in the region of small negative values for the next-nearest neighbour interac- 
tion. Monte Carlo simulations of the universal quantity of the Binder cumulant and 
the universal susceptibility ratio with different next-nearest neighbour interactions 
Jnnn showed a variation of both quantities for all applied models. The variation of 
these universal quantities leads to the conclusion, that in the considered region the 
NNN Ising model does not belong fully to the 3D Ising universality class. 
We assume that for regions, where chosen values of Jnn and Jnnn still cause fer- 
romagnetic phase transitions, but are close to a region where we will find no more 
ferromagnetic behaviour, a new definition for universality classes for this model has 
to be found. 
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